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We have investigated the reduction of the glass transition 
temperature, Tg, for thin supported films of particularly small 
molecular weight {Mw = 2 kg/mol) polystyrene, and found 
good agreement with earlier studies on larger molecules. By 
combining the eigenmodes of the films with mode coupling 
theory, we arrive at a simple and predictive model which is in 
quantitative agreement with the data. Its only fitting param- 
eter is the shear modulus governing the relevant fluctuations. 
The weak dependence of the latter upon the molecular weight 
of the polymer is in accordance with the general observation 
that the reduction of Tg is largely independent of chain length 
at sufficiently low molecular weight. 



The experimental observation that thin polymer films 
melt at temperatures well below the glass transition tem- 
perature of the bulk polymer Q is one of the major cur- 
rent challenges in the theory of homopolymers. Although 
one might at first glance expect such a behaviour, due 
to the impact of the finite size geometry of a thin film 
upon objects as large as polymer molecules, a slightly 
closer look reveals the intricacy of the problem: the ef- 
fect of a noticeable reduction of the glass temperature 
Tg, namely, occurs in many cases at film thicknesses or- 
ders of magnitude larger than the end-to-end distance, 
Ree, of the molecules. A fundamental understanding of 
this effect would be of great interest not only for polymer 
physics, but potentially also elucidate the physics of the 
glass transition in general. 

For the sake of clarity, but without loss of sufficient 
generality, let us restrict our discussion to the case of 
polystyrene (PS), which is among the best studied poly- 
mers in this context. As it has recently been pointed out 
in very thorough and elucidating studies [^|j3j, it is nec- 
essary to distinguish between long chain polymers with 
molecular weight above about 380 kg/mol, and shorter 
ones. For long PS chains, with a molecular weight of, 
say, 500 kg/mol, a substantial reduction is observed at a 
thickness below approximately 60 nm. Although Ree is 
still less than that, the length of a stretched molecule of 
this weight is of the order of one micron. In this case, 
the effect might indeed be attributed to finite size effects 
upon the chain conformation and chain motion. These 
ideas have been meanwhile developed and seem to repro- 
duce successfully the features experimentally observed 
i|. 

For shorter chains, however, the puzzle is more in- 



volved, since it is meanwhile accepted that the charac- 
teristic dependence Tg{h) upon the thickness is largely 
independent of the chain length |^,|| in this case. A 
finite size effect upon the geometry and motion of the 
single chain is thus ruled out as the basic mechanism. 
Many attempts have been undertaken to explain these 
findings on other grounds, mostly by considering micro- 
scopic models of the inner structure of the films. D. Long 
and F. Lequeux have envisaged the freezing of the film as 
a percolation of rigid domains |0] , mediated by thermal 
fluctuations Other models assume a layered struc- 
ture of the film, with a particularly mobile region close to 
the free surface of the film. Within this framework, For- 
rest and Mattsson 1^ have recently been able to achieve 
quite impressive accordance with the experimental data 
1^,^ . Their model makes use of the so-called cooperativ- 
ity length, ^(T), which plays a mayor role in a whole class 
of theoretical concepts of the glass transition. The only 
drawback is that there is yet no well established theory 
of ^{T). Furthermore, as a consequence of the two-layer 
structure of the film inherent in the model, it is not com- 
pletely clear why there should not be two glass transitions 
rather than a single one shifted in temperature. 

As a complementary approach, computer simulations 
of polymer films with free surfaces have recently been car- 
ried out, and most of the experimental findings were ob- 
tained qualitatively However, the polymer chains in 
these simulations were very much shorter than those used 
in the experiments (16 monomer units), well below the 
entanglement limit. It is therefore worthwile to study ex- 
perimentally the behaviour of polymers with short chain 
length, both to ease comparison with simulation and to 
explore the range of validity of the apparent indepen- 
dence of the shift of Tg on molecular weight. 

We have thus investigated the glass transition in thin 
supported films of PS with a molecular weight as small 
as 2 kg/mol (~ 20 monomer units). The polymer was 
purchased from Polymer Labs (UK) with a polydisper- 
sity index M\y /M^ = 1.05, the radius of gyration is 1.3 
nm. Effects from the molecular geometry are thus ex- 
pected only for films of few nanometers thickness. The 
films were spin cast from toluene solution onto silicon 
wafers (Silchem GmbH, Freiberg/Germany), which were 
previously cleaned by ultrasonication in acetone, ethanol, 
and toluene, subsequently. Residual organics were re- 
moved with Caro's acid (1:1 H2SO4 + H2O2), and the 
substrates were thoroughly rinsed with hot millipore wa- 
ter afterwards. Films were investigated with thicknesses 
ranging from 4 to 160 nm. The roughness of the free 
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surface of the films was less than 0.2 nm, as revealed by 
scanning force microscopy (SFM). 

The glass transition temperature was determined in 
several ways independently. The standard procedure of 
monitoring the thermal expansion of the film via ellip- 
sometry, as introduced by Keddy and Jones [Q , was used 
for film thicknesses down to 9.6 nm. Thermal expansion 
coefficients agreed well with previously published data 
both above and below Tg. Possible loss of material dur- 
ing the heating and cooling procedures was checked by 
applying the Clausius-Mosotti relation to the thickness 
and refractive index data obtained from the ellipsometry. 
Perfect conservation of material was invariably found. At 
large film thickness, Tg obtained in this way approached 
327 ± 1 K, which is consistent with the temperature at 
which macroscopic melting is observed in the bulk {Tg) 
for PS with this chain length. 

For thinner films, we observed the buildup of amplified 
thermal fluctuations (spinodal dewetting) |p^ , |lT| . These 
processes were monitored by SFM with in-situ heating. 
In order to speed up the experiments to a feasible time 
scale, dewetting was observed at temperatures close to 
Tg. In this way, the viscosity ry(T) was determined. Fur- 
thermore, we found that rj{T) varied with temperature 
according to a Vogel-Fulcher law, just as the bulk vis- 
cosity, but with the Vogel-Fulcher temperature shifted 
by a certain amount depending on film thickness. In this 
way, the apparent glass transition temperature of the film 
could be obtained by extrapolation. 

Our data are shown in fig. |l] as the full symbols. The 
circles represent the thermal expansion measurements, 
the squares were obtained by extrapolation from spin- 
odal dewetting experients. As one can clearly see, the 
glass transition temperature is substantially reduced for 
all films thinner than about 50 nm. The solid line repre- 
sents the function 



Tg = T"(l + h^/h)- 



(1) 



This form has been shown before to account well for the 
data obtained by others for larger molecular weight films, 
if /lo = 0.68 nm was assumed for PS ||l^. Within ex- 
perimental scattering, our data exhibit indeed the same 
dependence of Tg{h) in the full range of film thickness 
explored. Although our polymer chains are roughly by 
a factor of 50 shorter, we obtain /iq = 0.82 from the fit, 
which is quite close to the above value. This confirms 
the weak dependence (if there is a significant one at all) 
of the reduction of the glass transition temperature on 
the molecular weight of the polymer, down to a molecu- 
lar weight as small as 2 kg/mol, comparable to what has 
been studied by computer simulation Q . 

It is clear that the behaviour displayed in fig. |^ can in 
no way be attributed to the geometrical impact of the 
finite film thickness upon the microscopic conformation 
of the individual chains. Were therefore present here a 



novel approach to the problem which intentionally makes 
no reference to the molecular structure of the film. In- 
spired by the 'melting' of the spinodal waves, we con- 
sider instead the eigenmode spectrum of the (viscoelas- 
tic) polymer film. This can be considered with all possi- 
ble boundary conditions at the substrate, such that films 
with large extrapolation length, or free standing films, 
may as well be treated within the same framework. 

The spectrum can be obtained in a straightforward 
manner by combining standard theory of elasticity p^ ] 
and hydrodynamics p^ ] in the limit of small Reynolds 
number (Stokes dynamics). The equation of motion reads 



(2) 



where E is Young's modulus, 77 is the viscosity, and wo is 
the rate of relaxation of the individual chains by diffusion. 
is a vector field related to the strain tensor, S. If we 
restrict ourselves to one lateral {x) and one normal (z) 
coordinate, 6 is defined via 



S = 



dx4>x 



^{dx(l)z + dz(l>x) 
dz4>z 



(3) 



Finally, p is the pressure field. 

For harmonic excursions of the free surface, ({x) — 
Co exp{iqx — tot}, eq.(||) has solutions of the type 

02; = [1 + {h + <l~^)oi{<l)] cosh qz — q^^a{q) sinh qz 

0z = [1 + ha{q)] smhqz — za{q) coshgz (4) 

where for the function a(q), we find 



"(9) = (f ) ^ 



2qh _ 2 



1 + qh 



(5) 



for free standing films (symmetric modes) as well as for 
supported films with full slippage. At the free surface, 
we used the boundary conditions of zero tangential stress 
and p — —crdxxCi where a is the surface tension of the 
polymer. For free standing films, h is defined here as half 
the film thickness, such that the form of eqs. (|) and (|) 
accords with the experimental observation that the shift 
of Tg for a free film of thickness 2h is roughly the same 
as for a supported film of thickness h [ p^ . 
For the relaxation rates of the modes, we get 

2u; = {^0 + - + TT^) 
rj Zrjayq) 



, // E aq^ , . , 2ag2 
± J (^0 + - + TT^? + ^O^T 



(6) 



As it will become clear below, it is mostly the high fre- 
quency modes which contribute appreciably to the re- 
duction of the glass transition temperature. We thus 
consider only the high frequency branch of eq. (|^) . Since 
^0 « -^j the latter is approximately given by 
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w = cjQ H \- 



Tj 2r]a{q) 



(7) 



which is monotonously increasing with q. However, 
modes with q much larger than the inverse film thick- 
ness, , do not penetrate appreciably into the film, 
such that only a small fraction of the material takes part 
in these modes. Hence we are led to considering chiefly 
those modes for which q w /i"^, since these are the high- 
est frequency modes comprising more or less all of the 
film material. 

The central assumption is now that the physical cause 
for the melting, or freezing, of these thin film modes are 
memory effects in the polymer material. These are of 
course not included in the above linear theory, and may 
be formulated in a generic way by means of a suitable 
memory kernel, which may be written as m{0(i)} — 
ai<^-|-a2(/'^-f 030^-1-... It should be well noted that 

we are considering the strain field, 0, which describes the 
local state of the material. The equation of motion (2) is 
of first order in time, hence we have 

(P'+LU(j)+\[ m{(l>{T)}(l>'{t~T)dT (8) 

where A is a coupling parameter. This type of mode 
coupling equations has been thoroughly analysed [|8| in 
relation to the microscopic physics of the glass transition 
(to which we do not refer here), as well as to large scale 
degrees of freedom . It was found that the coefficients 
tti of the memory kernel vary concurrently with tempera- 
ture, and that upon crossing a certain border in the space 
spanned by the a^, the system freezes into a nonergodic 
state [0 . The existence of such a freezing transition has 
been found to be largely independent of the precise form 
of the memory kernel. For the simple scaling analysis to 
be done here, we therefore simply assume the Oi to vary 
all in the same way, which appears to be justified closely 
enough to the transition. In what follows, we thus as- 
sume that all temperature dependence may be absorbed 
into the coupling parameter. A, and the system freezes 
when a critical value, Ac, is reached. 

From eq. js] we see directly that A has the same dimen- 
sion as Lo, and there is no further independent time scale 
in the system. Since we expect A — > for T — > oo, and 
A ^ cxD for T — > 0, we can expand A = bi^+b2{-^^ + 
with E&i = Ac. Retaining only the leading term pO[ , and 
setting q — 1/h as discussed above, we arrive at a simple 
formula for the glass transition temperature: 



Tg(i?)=T" 1 



1.16 a 



(9) 



This is exactly the form of eq. (|l|), and fig. 1 shows that 
this describes our data very well. There was no significant 
improvement by including higher terms in the expansion 
series for A. 



The quantity rjcoo is on the order of a few kPa and can 
in general be safely neglected against E. The characteris- 
tic length scale which appears here, and which was called 
Hq in eq. (|l|), is thus the ratio a/E, with the surface ten- 
sion a — 31 mN/m for PS. It should thus be well noted 
that the elastic modulus determining the dynamics of the 
relevant modes, E, is the only physical fitting parameter 
in the model. 

The physical picture which emerges from this model 
is that as the temperature is increased, melting proceeds 
as the fastest mode involving all of the film material de- 
taches from its freezed state and fluctuates. The con- 
comitantly increased motion of polymer chains reduces 
the effective viscosity also for modes with smaller q, 
which, as a consequence, are sped up and melt in turn. 
In this way, the film finally melts at all length scales. One 
might apply this view also to the surface of a bulk poly- 
mer sample, for which the spectrum is obtained setting 
a = q/2. As it is readily seen, one should expect surface 
melting at the polymer surface, down to a thickness of 
hsm = hoT°{T^ - T)-'^. To the best of our knowledge, 
there is yet no conclusive evidence in favour or in dis- 
favour of polymer surface melting. This may be left to 
further studies. 

Let us turn back to our experimental data. From Hq = 
0.82 nm, as obtained from the fit, we find for the elastic 
modulus E — 44.1 MPa, which is, on logarithmic scale, 
right in between the modulus of the frozen material (a 
few GPa) and the modulus just above Tg (about 300 
kPa). As it is well known, both values display no marked 
dependence on the molecular weight, such that on the 
basis of our model, Tg{h) is expected as well to be largely 
independent of molecular weight. This is in agreement 
with all data collected so far for molecular weights up 
to about 300 kg/mol. Nevertheless, the exact physical 
significance of this somewhat arbitrary choice of E, which 
may be viewed as E atTg, is to be investigated in further 
studies. 
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FIG. 1. The glass transition of thin films of 2 kg/mol 
polystyrene, as determined from thermal expansion (circles) 
and from the growth of spinodal waves (squares). The solid 
curve represents our model, which has the elastic modulus 
governing the dominant modes as its only fitting parame- 
ter. It furthermore corresponds to what was found for larger 
molecular weight PS before |l^. Top: linear scale. Bottom: 
logarithmic scale, showing more details at small film thick- 
ness. 
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